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Abstract 

Interference phenomena are the source of some of the spectacular colors of animals and plants in 
nature. In some of these systems, the physical structure consists of an ordered array of layers with 
alternating high and low refractive indices. This periodicity leads to an optical band structure that 
is analogous to the electronic band structure encountered in semiconductor physics; namely, specific 
bands of wavelengths (the stop bands) are perfectly reflected. Here, we present a minimal model 
for optical band structure in a periodic multilayer and solve it using recursion relations. We present 
experimental data for various beetles, whose optical structure resembles the proposed model. The 
stop bands emerge in the limit of an infinite number of layers by finding the fixed point of the 
recursive relations. In order for these to converge, an infinitesimal amount of absorption needs to be 
present, reminiscent of the regularization procedures commonly used in physics calculations. Thus, 
using only the phenomenon of interference and the idea of recursion, we are able to elucidate the 
concepts of band structure and regularization in the context of experimentally observed phenomena, 
such as the high reflectance and the iridescent color appearance of structurally colored beetles. 



I. INTRODUCTION 



It is well known from the seminal works of Bloch that the energy levels of periodic 
structures such as semiconductors consist of a set of individual bands, separated by gaps 
where no energy levels are present^. This result has important consequences for the electronic 
properties of these materials, and is standard textbook material in any introduction to solid 
state physics and quantum mechanics. Remarkably, the colors of many butterflies, birds 
and fish rely on an analogous mechanism, the periodicity of a physical structure in this case, 
to achieve colour appearance arising from optical band structure. They exhibit optical stop 
bands, namely wavelengths in the UV-visible-IR range where electromagnetic waves have 
a high reflectance^-, leading to the striking colors of many of these animals. The same 
phenomenon is used in technology to make highly reflecting coatings, some of which have 
nearly 100 percent reflectance over a chosen range of wavelengths^. Bloch's derivation for 
the band structure, as elegant as it is, demands certain knowledge of wave theory, and is 
possibly not the most natural approach when calculating optical reflectance. It is the purpose 
of this work to supply an alternative derivation for calculating reflectance, applied here to a 
minimal, one-dimensional, periodic model, which invokes only the concepts of interference 
and recursion. This approach is intended to make it more accessible to advanced high-school 
and degree-level students by using recursion relations. 

The power of recursion relations was extensively described in the work of Ignatovich^, with 
various applications including the case of a periodic potential. Here, we use the recursion 
relations to show explicitly how an optical stop band is formed, using the symmetry of the 
structure to simplify the calculations. Although previous works by Rayleigh^ and Huxley^ 
describe methods for calculating the optical reflectance of a series of periodic layers we believe 
the method outlined here for the consideration of optical stop bands offers an alternative 
approach. While much effort has been made by these previous works to find an exact 
solution for a finite number of layers, here we take advantage of the convergence of the 
recursion relations to a fixed point in the limit of an infinite number of layers. Furthermore, 
we explain that if the calculation is made for a finite number of layers N, then a certain 
regularization is needed when taking the limit of iV — > oo, which has the physical meaning 
of an infinitesimal absorption. Such regularization is commonly used in many other, more 
advanced limit calculations in theoretical physics. The analysis we present serves multiple 
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purposes: it introduces the concept of optical band structure; it provides an explanation for 
the phenomenon of structural colors in biological and technological systems to complement 
existing models and it provides a good example where regularization is needed in a theoretical 
physics model. 



II. THE MODEL 



Figures [TJi to [Tb show photographs of three structurally coloured green beetles; To- 
rynorhina flammea chicheryi, Chrysochroa raja, and Gastrophysa viridula respectively. Re- 
flectance spectra collected close to normal incidence from their elytra, using a conventional 
experimental approach^, are presented in [Tg. The reflectance maxima shown in this graph, 
which are normalized relative to each other, arise due to the interference of light in the 
layered structures that are found close to the elytral surfaces of each beetle. Figures [TU toQJ 
show transmission electron micrographs of cross-sections through these beetles elytra; stain- 
ing during the sample preparation process^, gives imaging contrast between the different 
layers. The layers indices of refraction are in the region of 1.55 and 1.68^^. 

We now present a simplified model for the optical structure; an optical system that 
consists of iV structural unit cells, each consisting of two slabs of transparent material with 
two different indices of refraction. We will denote the thicknesses of these two slabs l\ and 
I2, and their refractive indices rii and n 2 . The system is illuminated by monochromatic light 
of wavelength A, coming from a medium of air (index of refraction n = 1). For simplicity, we 
analyze the case where the light propagation is perpendicular to the slabs. For this analysis 
we take the simple case of taking all the slab widths to be constant, and the light coherence 
length to be infinite, (i.e. there is interference, constructive or destructive, between waves 
which move arbitrarily far into the material). Despite an absorptive component in the 
material of biological multilayer structures^, we neglect the absorption in the structure 
for now. This is in some aspects a subtle point, however, which we will later revisit. The 
geometry is depicted in Fig. [21 



In a similar fashion to Refs. 



10 and 



111 ], we use the solution for N unit cells to find a 



solution for N + 1 unit cells. In the following sections we will derive these recursive relations, 
for completeness, and then proceed to find their fixed point. This gives rise to the optical 
band structure. 
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FIG. 1. Photographs of the beetles (a) Torynorhina flammea , (b)Chrysochroa raja and (c) Gas- 
trophysa viridula. TEM cross-sections of the multilayers responsible for these spectra are presented 
in (d)-(f) for T. flammea, C. raja and G. viridula respectively. (Color photographs are courtesy of; 
(a) Richard Bartz (b) Didier Descouens, (c) James Lindsey at Ecology of Commanster). Reflection 
spectra taken from the elytra of these three structurally colored green beetles, normalized with 
respect to each other, are shown in (g), (black line - T. flammea, green line - C. raja and red line 
- G. viridula). 
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FIG. 2. The one-dimensional, periodic optical structure that leads to the formation of optical stop 
bands. The dotted rectangle denotes the unit cell. 



A. Background material and notation 

It is known that the reflection coefficient of a wave incident on the interface between two 
materials with indices of refraction n\ and n 2 is given by 16 : r\^ 2 = nl+nt • Notice that 
when the wave comes from a lower refractive index material the relative phase changes by 
7T, since the coefficient is negative. For a wave with wavenumber k, the spatial dependence 
of the phase is given by e ±tkx for a wave propagating to the right/left. The transmission 
for the same scenario is given by: ti > 2 = 2wi . We shall denote the wavevectors in the two 
materials k\ and k 2 . 

B. Reflectance from a single layer 

Let us begin by considering the following basic structure, depicted in Fig. El 
It is a standard calculation to find the reflectance and transmittance of this structure, 
by considering the multiple reflections inside the slab n 2 , as is usually done in the context 
of analyzing the Fabry- Perot interferometer—. This leads to the following sums of complex 
numbers for the transmittance and reflectance: 

t = t^h^e^+^il + (r 2 ^e lk2h ) 2 + (r^e^ 2 ) 4 ] (1) 

r = ri _ 2 e ik ' h + t^r^We^^+^^fl + (r^e*^ 2 ) 2 + (r^e^ 2 ) 4 ]. (2) 
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FIG. 3. The unit cell of the structure of Fig. EJ 



Fig. H] shows an example of a 'diagram', corresponding to one of the possible paths of a 
wave in the optical structure, which together lead to the infinite sums described in Eqs. (JTJ) 
and (EJ). 



n 1 


n 2 


n 1 













h h h 

2 2 

FIG. 4. An example for the path of a wave in a single layer structure. In this case the wave 
makes 5 internal reflections before exiting the layer and being reflected back (a finite, small in- 
cident angle was used in the figure, to make the path clear, but throughout the paper we will 
assume the wave propagation is perpendicular to the surface). This corresponds to the term 
ii-s>2^2-5>ii2^ie^ fel ' 1+2fc2 ' 2 ' l (r2_ ! .ie* fc2 ' 2 ) 4 in the equation for the reflectance in Eq. ([2]). 

It is important to point out that in Eqs. ([I]) and fl2]) the transmittance coefficient, t, is 
the amplitude of the wave travelling to the right, of the form f e lk ( x ~ x r)^ w here x r = l\ + I2 
describes the position of the right-end of the single unit cell 'Fabry-Perot' interferometer. 
In many cases a different notation is used, where t is the prefactor of te lhx ; clearly t in the 
two notations will differ by a phase e lkXr . We will adhere to the above notation throughout 
this manuscript, as it will simplify the resulting equations. 
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Summing up the trigonometric series, we find that in the case of n = 1 unit cells,: 

ti^ 2 We i(fcl/l+M2) 



t 



+ 



1 - (r 2 -ae' 



(3) 
(4) 



1 - (r 2 ^ie ifc2 ' 2 ) 2 

If the transmission and reflection from a series of n unit cells is t n and r n then we can 
determine t n+ \ and r n+ i, the transmission and reflection coefficients for n + 1 such unit 
cells. This is done by effectively adding the extra slab in front of a black box of n cells, the 
optical properties of which we know. We know that a single unit cell has the transmission 
and reflection coefficients t and r, as given by Eqs ([3]) and (j3J). Now, we can consider the 
transmission and reflection through the whole structure, as illustrated in Fig. [51 Clearly, 
it has the same mathematical structure as the single unit cell Fabry-Perot setup, which we 
solved previously (it is even a simpler situation now, since in the absence of absorption, the 
reflectance and transmittance of the black box and the additional slab are clearly symmetric: 
this implies that there is no longer a difference between r 1 _ 5 . 2 and r 2 ^\ and between ti^ 2 
and £ 2 ->i)- 
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FIG. 5. Schematic depiction of the reflectance and transmittance properties of the system that 
lead to the recursion relations and which enable representation of its optical band structure. 

This leads to the following recursion equations, which are essentially a slight modification 
to the Fabry-Perot formulas of Eqs. (j3J) and PJ): 

t n +i — i (5) 

1 - rr n 

r t 2 

r n +i = r + — . (6) 

1 - rr„. 
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This recursive relation is at the heart of the derivation and enables exact calculation of 
the reflectance and transmittance of any finite number of unit cells or, as we shall shortly 
demonstrate, the limit of infinite unit cells. In that limit, we expect to reach a fixed point, 
where: 

tn+l =t n = tf (7) 

r„+i = r n = r f . (8) 

There are some subtle issues associated with this assumption on which we shall elaborate 
in section IIV1 In principle, Eqs. 03HS]) enable us to solve for the optical properties of an 
infinite array of unit cells. Although the equations are rather cumbersome, they can be 
significantly simplified for the case of a symmetric structure, as we choose to be the case 
here; every unit cell is symmetric by construction. In this case we can obtain a general 
relation between the reflectance and transmittance, which we shall now derive. Fig. [6^. 
illustrates the propagation of light through some symmetric structure, which has light waves 
with amplitude Iq = 1 incident upon it from both sides. The subsequent amplitudes A and 
B, which are complex, can be can be expressed in two ways: we can view A as the result of 
the reflection of the wave coming from the left, interfering with the transmitted wave from 
the right. This gives A = B = r + t. On the other hand, if we reverse time, we have to take 
the complex conjugate of all amplitudes, as described in Fig. [6b. Then we can view the 
amplitude 7=1 as the sum of the reflection of A* and the transmittance of B*. We therefore 
find: 



l = rA* + tB* = (r + t)(r* + t*). (9) 
Straightforward algebra gives us the useful relation: 










r 


2 




r 





So far, we have discussed the calculation of the optical properties of an extremely large 
number of unit cells; that is to say, the way in which light can propagate through them. 
Since this does not include the effect of light entering from outside such a black box, then we 
also have to consider the addition of a different cell (not identical to the unit cell). This will 
enable the calculation of the reflections in the interface between the surrounding medium 
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FIG. 6. Propagation of light through a symmetric structure (a), and its time-reversed path (b). 



of air and the periodic structure comprising the array of unit cells. This boundary region 
is schematically presented in Fig. [7J It is clear that in this case there are no multiple 
reflections, but this structure is not symmetric and so the reflectance and tranmittance are 
different when light is coming from the left (1) or from the right (r). It is evident that: 



(i - m)/(i + m), t^ r = 2e ik ^ 2 /(i + m), r r ^ 



(ni - 1)/(1 + ni), t P _+, 



2n 1 e lfcl ' 1/ ' 2 /(l +ni). Note the extra phase (exponential term) accumulated in traversing the 
distance l±/2 in a material with index of refraction n\. In order to connect this extra cell 
with the infinite structure whose optical properties we have already solved, we repeat the 
summation leading to Eqs. (j5J) and (EI), taking care to keep track of the asymmetry of the 
reflection and transmittance of this structure. This leads us to the following equations for 
the reflection r tota i and transmittance t tota i of the structure of Fig. [2] comprising a large 
number of periods: 



•'total 



f total — r l^>r + 



tl-H-tf 

1 - r r ^ir f 



1 — 7V 



>"7 



(11) 
(12) 
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FIG. 7. The final cell which has to be added to the left of the infinite array of unit cells shown in 
Fig. [3j The precise position of the right boundary of the cell is not important, since it goes in the 
middle of a region of the same index of refraction n\. However, the left boundary must be chosen 
left of the interface with the air, to account for the additional reflectance from it. For convenience, 
we choose it arbitrary close to it, to avoid having an additional length scale in the problem, which 
obviously will not affect the physics (but only add an additional phase to the overall reflectance). 



III. RESULTS 



We are now in a position to under stand the emergence of the optical stop bands, namely, 
the regimes of wavelengths for which there is 100 percent reflection. If the fixed point 
reflectance rf = 1, i.e., the infinite array has a 100 percent reflectance, it is clear that this 
must also be the case when we add the additional cell of Fig. describing the air interface: 
since waves cannot propagate through the infinite structure in this wavelength band and 
since there is no energy loss in the structure (i.e. the transmittance vanishes) they must all 
be perfectly reflected. Therefore, in order to understand the stop bands, it is sufficient to 
understand the case when rf = 1, without using Eqs. (ITT]) and (j!2p (which will be important 
if one is interested in the value of the reflection when it is not 100 percent). From Eq. (JBJ), 
we obtain: 



rd 2 

r f = r + -^ . (13) 

1 — rr/ 

This leads to a quadratic equation for r/. The condition for no transmission is |rj| = 1. 
Combining Eqs. ffTUj) and f lT3|) . one can simplify the quadratic equation to the form: 
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r) + { ^±p-r f + 1 = 0. (14) 

Therefore an explicit expression for the reflection of the structure (without taking into 
account the additional air interface) is: 




r / = -^i±j (r +;; )2 -i. (15) 

; 2r 2 \ 4r 4 



All the coefficients of Eq. (|14p are real, and it follows that |r/| = 1 if and only if the 
linear term is smaller in magnitude than 2. This reduces to the condition: 

Real[r] < \r\ 2 . (16) 

In this case both solutions of Eq. ( fl4l) give |r/| = 1. Otherwise, only one of the solutions 
has 1 77 1 < 1, which is therefore the correct solution. 

Taking r from Eq. @, the condition of Eq. ffl6|) can be simplified to: 

I cos(fciZi) cos(fc 2 Z 2 ) sin(fc 1 Z 1 ) sin(/c 2 / 2 )| > 1- (17) 

2nin 2 

The expression is symmetric with respect to n\ and n 2 , as it must be: if the array has 
zero transmittance, it cannot support the propagation of light waves through it, and this 
must be the case whether the first slab is rii or n 2 . Fig. [8] is a plot of r tota i versus wavelength 
(essentially a measure of fciA), obtained from Eqs. f fl5|) and f fT2|) . for the arbitrary choice of 
A; 2 i 2 = kill {i.e., l 2 = Ixni/n?) which Land^ describes as an ideal multilayer. Note that the 
derivation did not assume anything about the relative optical thickness of the two layers, 
and the choice of equal optical thickness used for Fig. [8] is arbitrary. We find that the 
reflectance is 100 percent over a specific band of wavelengths and not just at one specific 
wavelength - this corresponds to an optical band gap. It is analogous to an electronic band 
gap encountered in semiconductor science. 

For a typical biological example, taking realistic values of the real part of beetle cuticle 
materials' refractive indices, rii = 1.55, n 2 = 1.68, h = h = 0.09 microns, we can numerically 
find, from Eq. ( fl7|) . that the edges of the stop bands lie at wavelength of 0.567 microns and 
0.597 microns. Fig. [8] shows the reflectance as a function of wavelength, found from Eq. 
( !T4|) . where this stop band can be seen, as well as the non 100 percent reflectance outside 
the band. 
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FIG. 8. Result for the reflectance from an infinite array of unit cells, as found from Eqs. (|15p and 
(|12p . m = 1.55, n% = 1.68, I = I2 = 0.09 microns. The region of perfect reflectance is the optical 
stop band. The inset shows a plot of the reflectance for a broader range of wavelength (plotted in 
the inset in terms of the wavenumber), and demonstrates that there are in fact an infinite number 
of additional stop bands, outside the visible spectrum. 

The fact that the reflectance depends in a sensitive way on the wavelength, implies that 
the interference of waves through this structure (made only from two transparent materials) 
will give rise to color. This color depends on the wavelength and the system's structural 
parameters. In this way, such as system is said to produce structural colors, and is the 
underlying mechanism responsible for the colors of many insects, birds and fish. Moreover, 
since the reflectance is 100 percent over a range of wavelengths and the color reflected can 
be tuned or tailored through choice of materials and their spatial geometries, this type 
of structure can be used to make high reflectance coatings, which are necessary for many 
applications in optics. We should point out here that the minimal model which we have 
described and used will not on its own suffice to explain certain properties which are observed 
in real beetles and other animals and plants, such as the width of the reflection maxima 
and the value of the peak reflectance (which is smaller than unity). Among the reasons for 
this is that we have chosen in this introductory description, not to represent the genuine 
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optical absorption of realistic beetle material^ 1 ^ and that of other biological systems in the 
modeling and examples we present here. We maintain also that certain extents of relative 
disorder in such biological structures, imposed upon their overall periodic order, is crucial 
to quantitatively explain these current discrepancies. We plan to address these issues in 
subsequent works. 

IV. THE NECESSITY FOR REGULARIZATION 

By using the recursion relations (jSJ) and (Q, we can readily find numerically the reflectance 
from any finite number of layers. In fact, the equations can be easily modified such that at 
every stage of the recursion the number of layers is doubled (all one has to do is replace r 
and t by and in Eqs. (J5J) and (jBJ)). This means that using m iterations of the recursion 
relations we find the exact solution to the reflectance from 2 m layers. Fig. [9] shows the result 
of the reflectance from 32 layers obtained in this way. The result for the reflectance as a 
function of the wavevector shows clear oscillations, which are not reduced in amplitude as 
the number of unit cells increases. In other words, the recursion does not seem to converge 
to the solution that we discussed before! However, upon the addition of an infinitesimal 
amount of absorption (which implies using a wavevector with a small imaginary component, 
since the wave propagation is described by e lkx ), the oscillations are diminished, and for a 
large number of unit cells the reflectance asymptotically approaches the fixed point solution 
of Eq. (fl5|) . This is demonstrated in Fig. [TOj for 2 16 unit cells. 

This type of regularization is often needed in calculations in physics, such as the scatter- 
ing in the Born approximation from a Coulomb potential or the calculations of the Green 
function of various wave equations, where an imaginary part ie has to be added before taking 
the relevant limit. Also here, the order of limits is important: one has to add some absorp- 
tion ie (which is certainly present in reality), and take the limit iV — > oo, before taking the 
limit e — > 0. Using this order of limits, we recover the fixed-point solution of the recursion 
relations, associated with Eq. (I15j) . as shown in Fig. [TU1 

For a given, small, absorption, we expect to see a transition from the oscillating behavior 
of Fig. |9] to the fixed point solution, as shown in Fig. [101 This transition is shown in 
Fig. [HI where the reflectance as a function of wavelength was calculated for a structure 
with a number of layers varying between 2 and 2 16 , with a small absorption such that 
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FIG. 9. Result for the reflectance from 32 unit cells, for the same parameters as Fig. [8J The 
oscillations in reflectance on either side of the maximum reflectance position persist even as N — > oo. 
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FIG. 10. Result for the reflectance from 2 16 unit cells, with a small imaginary part added to the 
wavevector. The parameters are as in Fig. [HJ with a small imaginary component added to the 
wavevectors such that Im\k\li\ = Im^fa] = 10 . 

Imlfah] = Im[k 2 k) = 10" 4 . 

V. CONCLUSIONS 

We have shown how one can understand the formation of optical stop bands using only 
elementary tools, namely, an understanding of interference and the utilization of recursion 
relations. Upon calculating the reflectance from an infinite periodic optical structure, we 
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FIG. 11. The analysis of Fig. [10] is repeated for the same parameters, but with a number of layers 
varying between 1 and 1000. A transition is observed between oscillatory behavior for a small 
number of layers, where the small absorption used does not play a role, to a large number of layers 
where the absorption makes the solution converge to the fixed point solution, with the emergence 
of a stop band. 



have shown that there is 100 percent reflectance for certain ranges of wavelengths, a phe- 
nomenon that underpins the existence of optical band structure. Moreover, for finite periodic 
structures, we have shown that in order for the reflectance not to oscillate strongly outside of 
these stop bands, one must add a token small amount of optical absorption to the material's 
refractive index. 

Due consideration of optical band structure is important for understanding a variety of 
natural phenomena, such as the structural colors of many animals in nature, the construc- 
tion of high-reflectance coatings and photonic devices. It sits next to the band structure 
of semiconductor materials as a centrally important fundamental concept in fundamental 
modern physics. 
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